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Results of WEBT, VLBA and RXTE monitoring of 3C 279 during 

2006-2007* 

V. M. Larionov 1 ' 2 , S. G. Jorstad 1 ' 21 , A. P. Marscher 21 , C. M. Raiteri 3 , M. Villata 3 , 1. Agudo 4 , M. F. Aller 5 , 
A. A. Arkharov 2 , 1. M. Asfandiyarov 20 , U. Bach 6 , R. Bachev 7 , A. Berdyugin 28 , M. Bottcher 8 , C. S. Buemi 16 , 
P. Calcidese 9 , D. Carosati 10 , P. Chariot 11 , W.-P Chen 12 , A. Di Paola 13 , M. Dolci 14 , S. Dogru 15 , 
V. T. Doroshenko 40 ' 34 ' 41 , Yu. S. Efimov 34 , A. Erdem 15 , A. Frasca 16 , L. Fuhrmami 6 , P. Giommi 36 , L. Glowienka 17 , 

I A. C. Gupta 4418 , M. A. Gurwell 19 , V. A. Hagen-Thorn 1 , W.-S. Hsiao 12 , M. A. Ibrahimov 20 , B. Jordan 40 , M. Kamada 22 , 
Q ; T. S. Konstantinova 1 , E. N. Kopatskaya 1 , Y. Y Kovalev 6 ' 23 , Y A. Kovalev 23 , O. M. Kurtanidze 24 , A. Lahteenmaki 25 , 

L. Lanteri 3 , L. V. Larionova 1 , P. Leto 37 , P. Le Campion 11 , C.-U. Lee 26 , E. Lindfors 28 , E. Marilli 16 , 1. McHardy 27 , 
O ; M. G. Mingaliev 42 , S. V. Nazarov 34 , E. Nieppola 25 , K. Nilsson 28 , J. Ohlert 29 , M. Pasanen 28 , D. Porter 30 , T. Pursimo 31 , 

J. A. Ros 32 , K. Sadakane 22 , A. C. Sadun 33 , S. G. Sergeev 34 - 41 , N. Smith 39 , A. Strigachev 7 , N. Sumitomo 22 , 
^ ! L. O. Takalo 28 , K. Tanaka 22 , C. Trigilio 16 , G. Umana 16 , H. Ungerechts 43 , A. Volvach 35 , and W. Yuan 18 

6 

'_ ■ Context. The quasar 3C 279 is among the most extreme blazars in terms of luminosity and variability of flux at all wavebands. Its variations in flux 
' and polarization are quite complex and therefore require intensive monitoring observations at multiple wavebands to characterise and interpret the 
^3 observed changes. 

1 i , Aims. In this paper, we present radio-to-optical data taken by the WEBT, supplemented by our VLBA and RXTE observations, of 3C 279. Our 

goal is to use this extensive database to draw inferences regarding the physics of the relativistic jet. 
' Methods. We assemble multifrequency light curves with data from 30 ground-based observatories and the space-based instruments SWIFT 
^ ' (UVOT) and RXTE, along with linear polarization vs. time in the optical R band. In addition, we present a sequence of 22 images (with po- 
, larization vectors) at 43 GHz at resolution 0.15 milliarcsec, obtained with the VLBA. We analyse the light curves and polarization, as well as the 

spectral energy distributions at different epochs, corresponding to different brightness states. 

Results. We find that the IR-optical-UV continuum spectrum of the variable component corresponds to a power law with a constant slope of - 1 .6, 
xj ' while in the 2.4-10 keV X-ray band it varies in slope from -1.1 to -1.6. The steepest X-ray spectrum occurs at a flux minimum. During a decline 
f--^ , in flux from maximum in late 2006, the optical and 43 GHz core polarization vectors rotate by ~ 300°. 

I ■ Conclusions. The continuum spectrum agrees with steady injection of relativistic electrons with a power-law energy distribution of slope -3.2 
qq ' that is steepened to -4.2 at high energies by radiative losses. The X-ray emission at flux minimum comes most likely from a new component that 
^ ' starts in an upstream section of the jet where inverse Compton scattering of seed photons from outside the jet is important. The rotation of the 
• . polarization vector implies that the jet contains a helical magnetic field that extends ~ 20 pc past the 43 GHz core. 

. — i , Key words, galaxies: active - quasars: general - quasars: individual: 3C 279 

X ■ 
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- - ■ 1. Introduction with extremely rapid flux variability, by the EGRET detector 

of the Compton Gamma Ray Observatory and at longer wave- 

Blazars form a class of active galactic nuclei in which the spec- length s near a historical maximum in brightness dWehrle et all 

tral energy distribution is dominated by highly variable non- 1998) xhe y _ my flare was coincid ent with an X-ray outburst 

thermal emission from relativistic jets that point almost directly without any lag longer than x day HoweV er, while one can draw 

along the line of sight. The quasar-type blazar 3C 279, at redshift some important conclusions from such campaigns, the behaviour 

Z = 0.538 (Burbid ge & Rosenbergj [1965D, is one of the most in- of 3C 279 i s typically too complex to characterise and relate 

tensively studied objects of this class, owing to its pronounced conc i usive l y t0 physical aspects of the jet from short-term light 

variability of flux across the electromagnetic spectrum (by more curV es at a few wavelengths 
than 5'" in optical bands) and high optical polarizatio n (highest 

observed value of 45.5% in U band. lMead et al.| [l990). Because Ultra-high resolution observations with the Very Long 
of this, 3C 279 has been the target of many multiwavelength Baseline Array (VLBA) have demonstrated that EGRET- 
campaigns, mounted in an effort to learn about the physics of detected blazars possess the most highly relativistic jet s among 
the jet and the high-energy plasma that it contains. For exam- compact flat spectrum radio sources (Jorsta det al.L l2001a; 
pie, in early 1996 the source was observed at a high -y-ray state, Kelle rmann et al.Ll2004t) . This is inferred from the appearance of 
superluminally moving knots in the jet separating from a bright, 



Send offprint requests to: Y.M. Larionov, vlar@astro.spbu.ru compact, stationary "core" in the VLBA images. Apparent 

* The radio-to-optical data presented in this paper are stored in the speeds in 3C 279 range from Ac to 16c (e.g.,|Jorstad et aU|2004 

WEBT archive; for questions regarding their availability, please contact and references therein). From an analysis of both the apparent 

the WEBT president Massimo Villata (villata@oato . inaf . it). speeds and the time scales of flux decline of individual knots in 
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3C 279, iJorstad etail d2005l) derived a Doppler beaming factor 
6 = 24 + 6, a Lorentz factor of the jet flow T = 16 + 3, and an an- 
gle between the jet axis and line of sight « 2°. An analysis of 
the times of high y-ray flux and superluminal ejections of y-ray 
blazars indicates a st atistical connection between the two events 
dJorstadetaU 12001b). An intensive set of multi-waveband mon- 
itoring and VLBA observations of BL Lac confirms this con- 
nection between high-energy flares and superluminal knots in at 
least one blazar (iMarscher et al.Ll2008l) . 

Given the complexity of the time variability of nonthermal 
emission in blazars, a more complete observational dataset than 
customarily obtained is needed to maximize the range of con- 
clusions that can be dr awn concerning jet physics. To this end, 
IChatteriee et al.ld200 8) have analysed a decade-long dataset con- 
taining radio (14.5 GHz), single-colour (R-band) optical, and X- 
ray light curves as well as VLBA images at 43 GHz of 3C 279. 
They find strong correlations between the X-ray and optical 
fluxes, with time delays that change from X-ray leading opti- 
cal variations to vice-versa, with simultaneous variations in be- 
tween. Although the radio variations lag behind those in these 
wavebands by more than 100 days on average, changes in the 
jet direction on time scales of years are correlated with and 
actually lead long-term variations at X-ray energies. The cur- 
rent paper extends these observations to include multi-colour 
optical and near-infrared light curves, as well as linear polar- 
ization in the optical range and in the VLBA images, although 
over a more limited range in time (2006.0-2007.7). The radio-to- 
optical data presented in this paper have been acquired during a 
multifrequency campaign organized by the Whole Earth Blazar 
Telescope (WEBTflOur datase t includes the data fro m the 2006 
WEBT campaign presented by Bottch er et al. d2007l) as well as 
some of the data analysed by Chatteri ee et alj (|2008). The first 
of these papers found a spectral hardening during flares that ap- 
peared delayed with respect to a rising optical flux. The authors 
interpreted such behavi our in terms of inefficie nt particle accel- 
eration at optical bands. Chatterie e et al. (2008) came to a simi- 
lar conclusion based on the common occurrence of X-ray (from 
synchrotron self-Compton emission) variations leading optical 
synchrotron variations, with the latter involving higher-energy 
electrons than the former. 

Our paper is organized as follows: Section|2]outlines the pro- 
cedures that we used to process and analyse the data. Section [3] 
presents and analyses the multifrequency light curves, Sect. [4] 
discusses the kinematics and polarization of the radio jet as re- 
vealed by the VLBA images, and Sect. [5]gives the results of the 
optical polarimetry. In Sect. [6] we derive and discuss the time 
lags among variations at different frequencies, while in Sect. [7] 
we do the same for the radio-to-X-ray spectral energy distribu- 
tion constructed at four different flux states. In Sect.|8]we discuss 
the implications of our observational results with respect to the 
physics of the jet in 3C 279. 



2. Observations, data reduction and analysis 

Table [1] contains a list of observatories participating in the 
WEBT campaign, as well as the bands/frequencies at which we 
acquired the data for this study. Reduction of the data generally 
followed standard procedures. In this section, we summarise the 
observations and methodology. 



2.1. Optical and near-infrared observations 

We have collected optical and near-IR data from 21 telescopes, 
listed in Table Q] The data were obtained as instrumental mag- 
nitudes of the source and comparison stars in the same field. 
The finding chart can be found at WEBT campaign WEB-page. 
We used optical a nd NIR calibration of standard stars from 
iRaiteri et aD (Il998l) and lGonzalez-Perez. Kidger. & Martm-Luisl 
(120011) . 

We carefully assembled and "c leaned" the optical light 
curves (see, e.g., Villa ta et"aTL l2002h . When necessary, we ap- 
plied systematic corrections (mostly caused by effective wave- 
lengths differing from standard BVRcIc bandpasses) to the data 
obtained from some of the participating teams, to match the cal- 
ibration of the sources of data that use standard instrumentation 
and procedures. The resulting offsets do not generally exceed 
(POl-WCB in/? band. 



2.2. Radio observations 

We collected radio data from 8 radio telescopes, listed in TableQ] 
Data binning and cleaning was needed in some cases to reduce 
the noise. We report also VLBA observations performed roughly 
monthly with the VLBA at 43 GHz (,17mm) in both right and 
left circular polarizations throughout 2006-2007 (22 epochs). 
In general, all 10 antennas of the VLBA were used to reach 
ultra-high resolution, ~0.15 milliarsecond (mas), for imaging. 
However, at a few epochs 1-2 antennas failed due to weather 
or technical problems. We calib rated and reduced t he data in 
the same manner as described in Jorst ad et alj d2005i) using the 
AIPS and Difmap software packages. The electric vector posi- 
tion angle (EVPA) of linear polarization was obtained by com- 
parison between the Very Large Array (VLA) and VLBA inte- 
grated EVPAs for sources OJ 287, 1156+295, 3C 273, 3C 279, 
BL Lac, and 3C 454.3 that are common to both our sample and 
the VLA/VLBA polarization calibration lisQ, at epochs simulta- 
neous within 2-3 days (10 epochs out of 22). At the remaining 
epochs, we performed the E VPA calibration using the D-term 
method (Gom eTet all 12002). We modelled the calibrated im- 
ages by components with circular Gaussian brightness distribu- 
tions, a procedure that allows us to represent each image by a 
sequence of components described by flux density, FWHM an- 
gular diameter, position (distance and angle) relative to the core, 
and polarization parameters (degree and position angle). 



2.3. UV observations by Swift 

In 2007, the Swift satellite observed 3C 279 at 15 epochs from 
January 12 to July 14, covering bot h bright and faint s tates 
of the source. The UVOT instrument ( Roming et al.L 120051) ac- 
quired data in the V, B, U, UW1, UVM2, and UVW2 filters. 
These data were processed with the uvotmaghist task of the 
HEASOFT package remotely through the Runtask Hera facil- 
it)El with software and calibration files updated in May 2 008. 
Following the recommendations given by lPoole et all (120081) for 
UVOT photometry, we extracted the source counts from a cir- 
cle using a 5" aperture radius, and determined the background 
counts from a surrounding source-free annulus. Since the errors 
of individual data points are of the order of 0?1, we use daily 
binned data for further analysis. 



1 http://www.oato.inaf.i t/blazars/webt/ (se e 
Bottcher et al., 2005; Villata et al., 2006; Raiteri et al., 2007). 



e.g. 2 http://www.vla.nrao.edu/astro/calib/polar 
3 http : //heasarc . gsf c . nasa . gov/hera 



Larionov et al.: Monitoring of 3C 279 during 2006-2007 
Table 1. Ground-based observatories participating in this work. 2006 

May 



Observatory 


Tel. diam. 


Bands 




Radio 




1 12.5 - 


Crimean (RT-22), Ukraine 


22 m 


36 GHz 


~3. 


Mauna Kea (SMA), USA 


8x6rr£ 


230, 345 GHz 


U_ x 12.0 


Medicina, Italy 


32 m 


5, 8, 22 GHz 


o 


Metsahovi, Finland 


14 m 


37 GHz 


11.5 - 

13 - 

14 - 


Noto, Italy 


32 m 


43 GHz 


Pico Veleta, Spain 


30 m* 


86 GHz 


ffi 15 " 

16 - 


UMRAO, USA 


26 m 


5, 8, 14.5 GHz 


RATAN-600, Russia 


600 m 1 ' 


1,2,5,8, 11,22 GHz 


17 - 


Near- 


infrared 




18 - 
13 - 


Campo Imperatore, Italy 


110 cm 


J,H,K 


14 - 


Optical 




^ 15- 


Abastumani, Georgia 


70 cm 


R 


Armenzano, Italy 


40 cm 


B,V,R,I 


16 - 


Belogradchik, Bulgaria 


60 cm 


R 


17 - 


Bordeaux, France 


20 cm 


V 


13 - 


Catania, Italy 


91 cm 


U,B,V 


14 - 


COMU Ulupinar, Turkey 


30 cm 


V,R,I 


DC 

15 - 


Crimean (AP-7), Ukraine 


70 cm 


B,V,R,I 


16 - 


Crimean (ST-7), Ukraine 


70 cm 


B,V,R,I 


Crimean (ST-7, pol.), Ukraine 


70 cm 


R 


12- 


Lulin, Taiwan 


40 cm 


R 


13 - 


Michael Adrian, Germany 


120 cm 


R 


— 14 - 


Mt.Maidanak, Uzbekistan 


60 cm 


B,V,R,I 


15 - 


Osaka Kyoiku, Japan 


51 cm 


R 


16 - 


Roque (KVA), Spain 


35 cm 


R 


12 - 


Roque (LT), Spain 


200 cm 


R 


Roque (NOT), Spain 


256 cm 


U,B,V,R,I 


-3 13 - 


Sabadell, Spain 


50 cm 


R 


14 - 


Sobaeksan, South Korea 


61 cm 


B,V,R,I 


St. Petersburg, Russia 


40 cm 


B,V,R,I 


11 - 


Torino, Italy 


105 cm 


B,V,R,I 




Valle d' Aosta, Italy 


81 cm 


R,I 


X 12 - 


Yunnan, China 


102 cm 


R 




° Radio interferometer includinj 


I 8 antennas of diameter 6 m. 


13 - 
10 - 


Data obtained with the IRAM 30-meter telescope as a part of stan- 


dard AGN monitoring program. 






X. 11 - 



Ring telescope. 



2.4. X-ray observations 

We observed 3C 279 with the Rossi X-ray Timing Explorer 
(RXTE) roughly 3 times per week, with 1-2 ks exposure dur- 
ing each pointing. This is a portion of a long-term moni- 
toring program from 1996 until the present (Marscher, 2006; 
Chatterie e et alll2008l) . Fluxes and spectral indices over the en- 
ergy range from 2.4 keV to 10 keV were computed with the 
X-ray data analysis software FTOOLS and XSPEC along with 
the faint-source background model provided by the RXTE Guest 
Observer Facility. We assumed a power-law continuum spectrum 
for the source with negligible photoelectric absorption in this en- 
ergy range. 



3. X-ray, optical and near-IR light curves 

Figure Q] displays the X-ray, optical, and near-IR light curves 
of the 2006-2007 observing seasons. Analysis o f the 2006 data 
is giv en in a p revious WEBT cam paign paper (Bottcher et all 
120071) : see also lCollmar et all (120071) . 

During nearly the entire 2006 season, 3C 279 was at a mod- 
erately bright optical level of R * 14™5, with two noticeable 
downward excursions in January and April (MJD 53740 and 
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MJD - 50000 

Fig. 1. X-ray, optical and near-IR light curves of 3C 279 in 
the 2006-2007 observing seasons. Previously published data are 
marked with blue (light) symbols. 



53830). At the beginning of the 2007 season, 3C 279 was at a 
brighter level, R = 12'. n 8. Subsequently, over a 100-day period, 
the flux decreased to R — 16^2. Superposed on this declining 
trend, we observed a mild outburst around MJD 54150-54180, 
with amplitude « CT7. The minimum light level at MJD 54230, 
though no t a reco rd for 3C 279, was close to a level at which 
iPian et all d 19991) were able to detect the contribution of ac- 
cretion disc in UV part of 3C 279 spectral energy distribution. 
Unfortunately, we were not able to detect it due to lack of Swift 
UVOT data close to these dates. 

One of the campaign goals was to evaluate the characteristic 
time scales and amplitudes of intranight variability. We found 
that, both in high and low optical states, such rapid variabil- 
ity does not exceed (JVQ2 hour -1 , while night-to-night variations 
reached U! 1 1-0^2. This can be confirmed by visual inspection 
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cc 

CO 




Fig. 2. Colour vs. magnitude dependence of 3C 279; data refer- 
ring to 2006 and 2007 seasons are marked with different colours. 



E, 20 




10 15 
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of Fig[T] where the general features of the light curves are not 
masked by noise. 

The bluer-when-brighter behaviour is easily seen on Fig. [2] 
(although it is less evident in the 2006 data). The regression 
curve is B - R cc (0.05 + 0.01)7?. There are two plausible expla- 
nations for the origin of this colour variation: either we observe 
simultaneously a constant (or slowly varying) source, e.g., the 
host galaxy and accretion disc, along with a strongly variable 
source (in the jet) with a substantially different SED; or, alter- 
natively, the bluer-when-brighter trend is intrinsic to the jet, as 
found for other blazars, e specially BL Lac object s ( Villata et al.l 



2000 ; iRaiteri et all 120031 fo r 0716+714; IVillata et al.l 12004 
Papadakis rVillata. & Raiterill2007l for BL Lacertae). In the fol- 
lowing we consider the former supposition as better substanti- 
ated for 3C 279 behaviour during these observing seasons. 



3.1. Analysis of optical data 

Following the technique developed by Hagen-Thorn (see, e.g., 
Hage n-Thorn et alll2008l and references therein), let us suppose 
that the flux changes within some time interval are due to a sin- 
gle variable source. If the variability is caused only by its fluxQ 
variation but the relative SED remains unchanged, then in the 
«-dimensional flux space [F\ , F„) (n is the number of spectral 
bands used in multicolour observations) the observational points 
must lie on straight lines. The slopes of these lines are the flux ra- 
tios for different pairs of bands as determined by the SED. With 
some limitations, the opposite is also true: a linear relation be- 
tween observed fluxes at two different wavelengths during some 
period of flux variability implies that the slope (flux ratio) does 
not change. Such a relation for several bands would indicate that 
the relative SED of the variable source remains steady and can 
be derived from the slopes of the lines. 

We use magnitude-to-flux calibration consta nts for op- 
tical (BVRI) and NIR (JHK) bands from iMead et alj 
(fl990lh and for U VOT (V, B, U, Wl, Ml, W2) - from 
Pool e et al.l (I2008"). The Galactic absorption in the direction 
of 3C 279 was calculated according to Cardelli's extinction 







B . 






U 






W1 












W2 



8 10 12 14 16 



F v [mJy] 

Fig. 3. Flux-flux dependences for optical and NIR (top panel) 
and UVOT bands (bottom panel). 



4 For the sake of brevity, we use the term "flux" instead of the more 
proper "flux density." 



law (ICardelli. Clavton. & MathisL Il989i) and A v = (T095 
dSchlggel, Finkbein er"& Davislll998T) ~ 

We use the better-sampled /?-band data to obtain relations in 
the form of Fi vs. Fr. Similarly, in the case of the UVOT data, 
we derive Fi vs. Fy. 

Figure [3] (upper panel) shows that the method described 
above holds true: during the entire time covered by our observa- 
tions, the flux ratios follow linear dependences, F , = A, + fi,- • Fr. 
Values of fi,-, along with the slopes of the regressions, can be 
used to construct the relative SED of the variable source. 

A similar analysis can be carried out with UVOT data. 
Figure [3] (bottom panel) demonstrates that the same kind of be- 
haviour is seen in the ultraviolet part of the 3C 279 spectrum. 
This confirms our supposition about self-similar changes in the 
SED of the variable source. 

Having found the slopes of flux dependences in the near-IR, 
optical, and UVOT ranges, we are able to construct the relative 
SED of the variable source in 3C 279. These results are given 
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Table 2. Multicolour properties of the variable optical emission 
component. 



Band 


log v (Hz) 




N 


r 


logfi 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Ground-based observations 


K 


14.1307 


0.010 


28 


0.989 


0.8286 


H 


14.2544 


0.016 


25 


0.998 


0.6688 


J 


14.3733 


0.026 


10 


0.998 


0.4609 


1 


14.5740 


0.055 


169 


0.998 


0.1520 


R 


14.6642 


0.076 






0.0000 


V 


14.7447 


0.095 


185 


0.998 


-0.1180 


B 


14.8336 


0.123 


121 


0.996 


-0.2434 


UVOT observations 


V 


14.7445 


0.093 






-0.1180 


B 


14.8407 


0.123 


14 


0.995 


-0.2607 


U 


14.9329 


0.147 


14 


0.999 


-0.4045 


Wl 


15.0565 


0.195 


13 


0.999 


-0.6438 


M2 


15.1286 


0.285 


6 


0.949 


-0.7306 


W2 


15.1696 


0.271 


14 


0.998 


-0.7856 



Columns are as follows: (1) - filter of photometry; (2) - logarithm 
of the effective frequency of the filter; (3) - extinction at the corre- 
sponding wavelength; (4) - number of observations in a given band, 
(quasi)simultaneous with R (for ground-based observations) and V (for 
UVOT observations); (5) - correlation coefficient; (6) - logarithm of the 
contribution of the variable component in the SED at frequency v rela- 
tive to the contribution at R band. 



1.0 




-1.0 1 1 1 1 1 1 1 1 

14.0 14.2 14.4 14.6 14.8 15.0 15.2 15.4 

log V [Hz] 

Fig. 4. NIR-optical-UV relative spectral energy distribution of 
3C 279 variable source(s), normalized to R band. 



in Table [2] (slopes in the UVOT bands are obtained relative to 
the UVOT V band and fitted to ground-based data) and shown in 
Fig- H] (The double B-band point is caused by different effective 
wavelengths for the ground-based telescopes and that of UVOT.) 
We note that we are unable to judge whether there is only one 
variable source acting throughout our observations or a number 
of variable components with the same SED. In any case, in this 
wavelength range we derive a power-law slope with a — 1.58 + 
0.01, in the sense F v oc v~ a . 

It is possible to "reconstruct" the light curves in the BVIJHK 
bands using the dependences shown in Fig. [3] This allows us to 
( 1 ) check the reliability of this procedure and (2) form an impres- 
sion of the variability behaviour without gaps that otherwise are 
present in less well-sampled bands. Figure [5] clearly shows that 
the actual data (red circles) closely match the calculated values 
(black dots). 





13 
14 










CO 


15 
16 
17 
18 


• 1 1 1 1 


I I 


, A 

mi 


Hi" i i i i r 




13 










> 


14 
15 
16 
17 


V 




v A 


S„. ■ ' 




12 






mi 


1 H-t-t-H 




13 
14 
15 
16 






, A 


M 4 :' 

i i i i i m \ \ \ i i i r 




11 
12 
13 







i M i i i 

.-' V 






14 












10 






M|M 

\ 

r v 


1 1 1 M | 1 1 | 1 1 1 1 




11 
12 
13 












8 
9 
10 
11 
12 
13 






i i i i 

■X 

- y 


1 INI 



3700 3800 3900 4000 4100 4200 4300 



MJD - 50000 

Fig. 5. Comparison of observational data (red circles) and cal- 
culated values according to dependences shown in Fig. [3] (black 
dots). 



Additionally, Figs.[3]and|5]allow us to conclude that there is 
no noticeable lag between any of the optical and near-infrared 
light curves; otherwise, we would see both hysteresis loops in 
Fig. [3] and misfitting in Fig. [5] If any lags in that wavelength 
range do exist, they are s horter than a few hou rs. This contradicts 
the inference drawn by Bot tcher et al.l d2007l) based on a more 
limited dataset. 

3.2. X-ray vs. optical fluctuations 

From an analysis of the p ower spectral density of longer- 
term monitoring observations, Chatteri ee et al.l (|2008) have con- 
cluded that the X-ray light curve of 3C 279 has a red-noise na- 
ture, with fluctuations on short time scales having much lower 
amplitude than those on longer time scales. In general, the X-ray 
variations ar e either less pronounced or similar to those at the op- 
tical R band. lChatteriee et al.l d2008l) explain this as the result of 
the higher energy of electrons that emit optical synchrotron ra- 
diation compared to the wide range of energies - mostly lower - 
of electrons that scatter IR and optical seed photons to the X-ray 
band. However, we see the opposite occur between MJD 53770 
and 53890, when the X-ray fluctuations were more dramatic than 
those at any of the optical bands (see Fig . [TJ . (Note: the apparent 
rapid X-ray fluctuations during the minimum at MJD 54210 to 
54240 is mostly noise, since the logarithm of the uncertainty in 
the fluxes during this time was of order +0.07.) We can rule out 
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an instrumental effect as the cause of the rapid X-ray variabil- 
ity, since the X-ray flux of PKS 15 10-089, measured under very 
similar circumstances, displays only more modest fluctuations. 
Although a transient, bright, highly variable X-ray source within 
~ 20' of 3C 279 could also have caused the observed fluctua- 
tions, no such source has been reported previously. Hence, we 
accept the much more likely scenario that the variations are in- 
trinsic to 3C 279. We discuss a possible physical cause of the 
rapid X-ray variations in Sect. 18.11 

4. Kinematics and polarization of the radio jet 

The Boston University group observes the quasar 3C 279 with 
the VLBA at 43 GHz monthly (if dynamic s cheduling works 
prope rly) in a program that started in 2001 (see lChatteriee et all 
2008). Figure[5]shows the sequence of total and polarized inten- 
sity images of the parsec-scale jet obtained with the VLBA dur- 
ing 2006-2007. The da ta were reduced and modelled in the same 
manner as described in Jorsta d et alj d2005l) . The images are con- 
volved with the same beam, 0.38x0.14 mas at PA=-9°, corre- 
sponding to the average beam for uniform weighting over epochs 
when all 10 VLBA antennas were in operation (16 epochs out 
of 22 epochs shown). The images reveal motion of two new 
components, C23 and C24 (we follow the sc heme of component 
designation adopted in lChatteriee et alj 2008). Figure 7 plots an 
angular separation of the components from the core vs. epoch. 
Although there is some deviation from ballistic motion (espe- 
cially for C23 within 0.3 mas of the core), a linea r dependence 
fits th e data according to the criteria adopted by Jorsta d~et alj 
(120051) . This gives a high apparent speed for both components, 
16.5 + 2.3 c and 14.7 ± 0.9 c, respectively, for C23 and C24 (for 
cosmological parameters H a = 70 km s _1 Mpc -1 , D. m - 0.3, 
£Xi = 0.7), and yields the following times of ejection of compo- 
nents (component's coincidence with the core): MJD 53888+55 
and 54063+40. In projection on the plane of the sky, these 
components move along similar position angles, PA ~ -114° 
and PA ~ -120°. The directions are different from the posi- 
tion angles of compone nts ejected in 2003-2004, PA ~ -150° 
dChatteriee et all [2008). when both the optical and X-ray activ- 
ity of the quasar was very modest. The images in Fig. [6] con- 
tain another moving component about 1 mas from the core (C21, 
Chatteriee et al. 2008), which is a "relic" of the southern jet di- 
rection seen near the core during 2004-2005. Using the method 
suggested in iJorstad et alj ([2005), we have estimated the vari- 
ability Doppler factor, 5 var , from the light curves and angular 
sizes of C23 and C24, 5 var = 25 and 5 var = 29 respectively, with 
uncertainties ~ 20%. From these values of 5 var and /3 app , we de- 
rive the Lorentz factor and viewing angle of the jet, F = 18 + 5, 
= 2.1°±O.3°forC23,andF = 18 + 3,0= 1.6° ± 0.4° for C24 
These values are consistent with those derived by Jorsta cfet alj 
d2005l) during the period 1998-2001 when 3C 279 was in a very 
active X-ray and optical state, similar to that seen in 2006-07. 

The polarization behaviour in the VLBI core at 7 mm, com- 
bined with the jet kinematics discussed above, helps to define the 
relationship between the physical state of the jet and the multi- 
waveband variations that we have observed. Figure [8] shows the 
light curve of the VLBI core at 7 mm and also degree and po- 
sition angle of polarization in the core. The dotted vertical lines 
indicate the times of ejection of components C23 and C24. The 
ejections of both components in 2006 occurred near maxima in 
the light c urve of the VLBI core re gion, in keeping with the con- 
clusion of lSavolainen et all ([2002 ) that the emergence of new su- 
perluminal knots is responsible for radio flares observed at ~ 40 
GHz. Each of the ejections in 2006 coincides with the start of a 
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Fig. 6. Total intensity VLBA images of 3C 279 at 43 GHz. Line 
segments within each image indicate the direction of polariza- 
tion electric vectors, with length proportional to the polarized 
intensity. The resolution beam is shown by the cross-hatched el- 
lipse in the lower left corner. The contour levels correspond to 
0.25,0.5,1,2,4,8,16,32, and 64% of the peak intensity of 17.7 Jy 
beam 1 . The solid lines indicate the positions of components 
C23 and C24 that were ejected in 2006. 



rotation of the EVPA in the core of duration ~60-80 days. During 
the rotation the degree of polarization in the core reaches a min- 
imum (< 0.5%). At the end of the rotation, the EVPA in the core 
aligns with the jet direction (see also Fig. [9]). 

Figure [8] also presents the X-ray and radio (37 GHz) light 
curves. The ejections of the components occurred just after max- 
ima in the X-ray light curve. (Although the peak of the second 
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Fig. 8. X-ray and radio observations of 3C 279 during 2006- 
2007. Solid lines mark VLBI core radiation at 7 mm, triangles - 
X-ray flux (arbitrary units) and diamonds - 37 GHz flux. Vertical 
dotted lines indicate the times of ejection of components C23 
and C24 seen on the VLBA images, and horizontal bars in the 
bottom panel - uncertainty of the times' determination. 



flare is not covered by our observations due to the proximity of 
the quasar to the Sun, we can infer the approximate time when 
it occurred by interpolation of the rising and declining phases of 
the event.) The X-ray maxima lead the ejections by ~60-80 days. 
The similarity of the behaviour at different wavelengths strongly 
suggests that the two major flares observed from X-ray to radio 
wavelengths are associated with disturbances propagating along 
the jet. 

5. Optical Polarimetry 

All optical polarimetric data reported here are from the 70cm 
telescope in Crimea and the 40cm telescope in St. Petersburg, 
both equipped with nearly identical imaging photometers- 
polarimeters from St. Petersburg State University. Polarimetric 
observations were performed using two Savart plates rotated by 
45°relative to each another. By swapping the plates, the observer 
can obtain the relative Stokes q and u parameters from the two 
split images of each source in the field. Instrumental polarization 
was found via stars located near the object under the assumption 
that their radiation is unpolarized. This is indicated also by the 
high Galactic latitude (57°) and t he low level of extinction in the 
direct ion of 3C 279 (A v = 0?095; ISchlegel. Finkbeiner. & Davisl 
I1998I) . 

The results of polarimetric monitoring of 3C 279 during the 
2007 campaign are given in Fig. [9] with /?-band photometry 
shown in the upper panel for comparison. The most remark- 
able feature of the polarimetric behaviour is the smooth rota- 
tion of the electric-vector position angle (EVPA) of polarization 
9r from the start of the observing season. This rotation contin- 
ued, while gradually slowing, over approximately two months. 
In total, the EVPA rotated by ~ 300°, eventually stabilizing at 
G R -240°. Note that, because of the +180°n ambiguity in G R , this 
value is the same as -120°, the direction of the VLBA inner jet 
in 3C 279 during 2007 (cf. Fig.©. The bottom panel of Fig. g] 
also shows the 7mm radio EVPA (red open circles), which ro- 
tates along with the optical EVPA, although with fewer points 
to define the smoothness of the rotation. The degree of optical 
polarization pu varies more erratically than the EVPA, reaching 
minimum values of ~ 2% close to the mid-point of rotation of 6r 
and varying between 5% and 30% during the remainder of the 
observational period. 

6. Radio light curves and cross-correlation analysis 

The common way to estimate delays between light curves at 
different wavelengths is to use the DCF (Discrete Correlatio n 
Function; lEdelson & Krolikll 19881; iHufnagel & Bregman|[l992l) . 
In the case of inter-optical delays, we have mentioned above 
that there are no lags between any pairs of wavelengths (see 
Sect. 13. p . This conclusion is confirmed with DCF analysis us- 
ing different pairs of optical bands - in all the cases the best fit 
is obtained for days lag. We failed to evaluate DCF lags be- 
tween optical and radio frequencies for both observational sea- 
sons simultaneously, because of the seasonal gap in optical ob- 
servations caused by solar conjunction and, to a lesser extent, the 
different time sampling. However, inspection by eye of the op- 
tical R and radio light curves (Fig. ITOb reveals progressive lags 
that increase toward lower frequencies, with a time delay be- 
tween the minima at R band and 37 GHz in mid-2007 of ~ 60 
days. The dashed slanted lines in Fig. [TUl connect the positions 
of two minima and one maximum that are easily distinguished 
in the light curves. The question mark in the /?-band panel de- 
notes the position of a tentative optical maximum, missed due to 
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Fig. 9. Optical polarimetric data for 3C 279 obtained during the 
2007 campaign. EVPA data for the 43 GHz core are marked with 
red (open) symbols in the bottom panel. The +180°« ambiguity 
in the EVPAs is resolved by requiring the change across contigu- 
ous epochs to be < 90°. No Faraday rotation correction has been 
applied to the 43 GHz data; see the text. 
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the seasonal gap. Such behaviour is not unexpect ed and confirms 
the results obtained durin g earlier epochs (e.g jTornikoski et all 
ll994[lLindfors et al.ll2006l) . 

We calculated the optical-radio (37 GHz) delays separately 
for two observing seasons (Fig. [TTb). The negative values of 
delays agree with our conclusion made from eye inspection 
of Fig. [10] In order to evaluate time lags quantitatively at ra- 
dio frequencies, we have computed the DCF for the 22 GHz, 
14 GHz, 8 GHz and 5 GHz light curves relative to 37 GHz. 
The results are plotted in Fig. [TTb , where vertical bars indi- 
cate the lags corresponding to the DCF centroids. Despite the 
large uncertainties, this plot confirms our conclusion about vari- 
ations at lower frequencies lagging be hind those a t high fre- 
quencies. Simi lar results we r e found bv lRaiteri et aL[ (2003) for 
S5 07 16+714, IVillata et all d2004bl) for BL Lac, iRaiteri et al l 
J2005h for AO 0235+ 164. lviilata et all J2007h for 3C 454.3. 

The upper panel of Fig.[TT]displays the X-ray-optical DCF. 
The correlation is rather weak, with changes in the optical (R) 
flux lagging behind X-ray variations by ~ 1 day. Inspection by 
eye suggests that the correlation is strong during 2007 and weak 
or non-existent in 2006. The most salient feature in the light 
curves, a deep minimum between MJD 54220 and 54250, oc- 
curs essentially simultaneo usly at both X-ray and optical wave- 
bands. As is discussed by Chatt eriee et al.l d2008l) . the X-ray- 
optical correlation fluctuated slowly over the years from 1996 to 
2007, both in strength and in the length and sense of the time 
lag, with the weakest correlation occurring when the lag is near 
zero. The latter behaviour matches that observed in 2006-07. 
The controlling factor may be slight changes in the jet direction 
over time, which modulates the Doppler factor, light-travel de- 



Fig. 10. Radio and 7?-band light curves of 3C 279; dashed slanted 
lines trace mimima and maxima of light curves at different fre- 
quencies. Error bars are shown only for 1 and 2 GHz data, where 
they are larger than symbol size. 



lay across the source, and the position in the jet whe re the main 
observed emission originates dChatteriee et al 1 120081 



7. Spectral energy distribution 

Figure Q~2] displays the spectral energy distribution (SED) of 
3C 279 at four different flux levels at epochs when we have suf- 
ficient data. Note that the non-variable optical component has 
not been subtracted in these plots. The overall increase in flux 
across the entire sampled frequency range strongly suggests that 
the nonthermal emission at radio through X-ray wavelengths has 
a common origin, which we identify with the jet based on the 
appearance of new superluminal knots associated with each of 
the two major outbursts in flux. The other striking feature of 
the SEDs is the steepening of the X-ray spectrum at lower flux 
levels. (Note: the X-ray spectral index is derived from an expo- 
nential smoothing of the values obtained from the observations, 
with a smoothing time of 10 days. This is necessary because 
of the large uncertainties in individual spectral indices derived 
when the flux is low.) It is perhaps significant that the highest 
X-ray spectral index of 1.6 matches that of the variable optical 
component, while the flattest X-ray spectrum has an index that 
is lower by 0.5. As we discuss below, this has important implica- 
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Fig. 11. Discrete correlation functions (a) for optical-X-ray, (b) 
optical-radio (37GHz) and (c) radio frequencies relative to the 
37 GHz light curve. Negative time lag corresponds to lower- 
frequency radiation changing later than at higher frequencies. 
Vertical bars in panel (c) indicate the lags corresponding to the 
DCF centroids. Note that in panel (b), due to seasonal gap in 
optical range, the DCF is calculated separately for the 2006 and 
2007 seasons. 



tions for the energy gains and losses experienced by the radiating 
electrons. 



8. Discussion 

8.1. Variability in flux and continuum spectrum 

During the two observing seasons reported here, 3C 279 ex- 
hibited flux variations from radio to X-ray wavelengths (see 
Fig. [TJ. We observed the shortest time scales of variability, 
At I \n{F m ax/f min) ~ days, in the X-ray and optical bands. There 
were two deep minima in the X-ray light curve, the first (at the 
start of 2006) coinciding with a low optical state with at least 
two closely spaced minima, and the second essentially simulta- 
neous with a deep optical minimum between MJD 54210 and 
54240. The discrete cross-correlation function indicates that the 
time lag between the X-ray and optical bands is essentially zero. 

The SEDs in Fig. Q~2] demonstrate that the X-ray emission is 
not simply a continuation of the o ptical synchrotron spec trum. 
This agrees with the conclusion of Chatter iee et al.l d2008l) that 
the X-rays result from synchrotron self-Compton (SSC) scatter- 
ing. During the highest flux state, the X-ray spectral index was 
a x = 1.1, while the optical spectral index of the variable com- 
ponent was Copt = 1 .6 throughout both observing seasons (see 
Fig. HI). This can be accommodated within a standard nonther- 



13.5 



13.0 



12.5 



10.5 



JD 2454114 
JD 2454233 
JD 2454266 
JD 2454290 



V 



\ 



v 



10 12 14 16 

logv [Hz] 



18 



20 



Fig. 12. Radio-to-X-ray spectral energy distributions at various 
epochs in 2007. For each epoch, a cubic spline is drawn through 
the points. The errors of slope in the X-ray part of SED are of 
the order of 0. 1 . 



mal source model if there is continuous injection of relativistic 
electrons with a power-law energy distribution of slope -3.2 that 
is steepened to -4.2 at high energies by radiative energy losses. 
This implies that there is a flatter spectrum at mid- and far-IR 
wavelengths, below ~ 10 14 Hz, with qtr = 1.1 that matches 
that of the high flux state. At minimum flux, the X-ray spectrum 
steepens to a x = 1.6 = a opt . 

The steep X-ray spectrum during the flux minimum of 2007 
implies that the entire synchrotron spectrum above the spectral 
turnover frequency had a steep spectral index of 1.6. In the con- 
text of the standard nonthermal source model mentioned above, 
essentially all of the electrons suffered significant radiative en- 
ergy losses at the time of the low state. That is, their radiative 
loss time scale was shorter than the time required for them to 
escape from the emission region. For this to be the case, the in- 
jected energy distribution must have a low-energy cutoff that is 
far above the rest-mass energy. We can match the time scale of 
optical and X-ray variability of 3C 279 if we adopt a magnetic 
field strength in the emission region of B ~ 1 G. Electrons ra- 
diating at the UV wavelengths must then have a Lorentz fac- 
tor y = (E/mc2) ~ 6000B" I/2 (5/25)" I/2 , those radiating at a 
break frequency of ~ 10 14 Hz have y ~ 1500B- I/2 (<5/25r 1/2 , 
and those emitting at 10 12 Hz have y ~ 150B" 1/2 (5/25r 1/2 . 
Given the turnover in the SED at ~ 10 12 Hz (cf. Fig. [T2l . we 
assume that either the variable emitting component becomes 
self-absorbed near this frequency or that there is a lower-energy 
cutoff (or sharp break) in the injected electron energy distribu- 
tion at y ~ 150B _1 ' 2 (<5/25) -1 ^ 2 . The time scale for energy losses 
from synchrotron radiation, as measured in our reference frame, 
is ~ 4fi- 3/2 (c5/25r 1/2 days at ~ 10 12 Hz if we adopt a Doppler 
factor of 25 (see Sect.[4]above). This is similar to the minimum 
X-ray variability time scale that we observed. The time scale for 
synchrotron energy losses allows the variability at optical bands 
to be as fast as ~ 3 hours as long as the emission region is smaller 
than ~ 5 x 10 15 (5/25) cm. 

The 2.4-10 keV X-ray emission results from a range of 
seed photons scattered by electrons having a range of energies 
(McHard v et al.L 119991) . If there is a break in the spectrum, as 
we infer, then the bulk of the X-ray emission is from seed pho- 
tons at frequencies below the break - v ~ 10 12 — 10 14 Hz - scat- 
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tered by electrons with energies below the break - y ~ 150— 
1500. The rapid optical variability in the low flux state implies 
that the magnetic field remained of similar strength as for the 
high state. However, the break frequency must have decreased 
to nearly 10 12 Hz for the SSC X-ray spectrum to have steep- 
ened to a x = a op t = 1.6. This could occur simply from adi- 
abatic cooling if the emission region expanded by a factor of 
10 between the high and low states. However, in this case the 
magnetic field would have dropped, increasing the time scale of 
variability and decreasi ng the flux by more tha n the observed 
factor of ~ 10 (cf., e.g., Mar scher & Gearl [l985l) . Alternatively, 
the time for electrons to escape the emission region could have 
in creased as the flux declin ed. For example, in the shock model 
of iMarscher & Gearl d 1985b . the electrons are accelerated at the 
shock front and advect toward the rear of the shocked region un- 
til they encounter a rarefaction at a distance x behind the shock. 
If the primary energy loss is due to synchrotron radiation, the 
break frequency v* °c x 2 B' 3 . A factor of 10 increase in the ex- 
tent (in the direction parallel to the jet axis) of the shocked region 
would then cause the break frequency to decrease from ~ 10 14 to 
10 12 Hz. However, this would require a high relative velocity be- 
tween the shock front and the rear of the shocked region, imply- 
ing that the emission near the latter is relatively poorly beamed 
and therefore not prominent. 

Another possibility, which seems less contrived, is that the 
X-ray emission during the flux minimum is not the decayed 
remnant of the previous event, but rather comes from the new 
component that caused the outburst seen during the last ~ 100 
days of the 2007 observing season. This component would first 
be seen well upstream of the location where the flux reaches 
a maximum, and therefore the external photon field from the 
broad emission-line region an d dusty torus would initiall y be 
higher than at later times (see ISokolov & M arscher. 2005|). In 
this case, the break frequency would increase with time as the 
component moves downstream and the level of inverse Compton 
energy losses decreases. This scenario predicts that there should 
be a y-ray flare from inverse Compton scattering of seed photons 
originating outside the jet prior to the main outburst. Although 
the y-ray observations needed to test this were not undertaken 
in early May 2007 (the newly-launched AGILE telescope was 
in the testing phase), future observations involving GLAST and 
AGILE will determine whether such precursor flares actually oc- 
cur. 

The higher amplitude of variability at X-ray energies relative 
to the optical bands between MJD 53770 and 53890 could be 
the result of the higher sensitivity of the SSC vs. the synchrotron 
flux to changes in the number of relativistic electrons (quadratic 
vs. linear dependence). Although, as discussed above, the X-rays 
during high flux states are scattered seed photons that originally 
had a range of frequencies from ~ 10 12 to 10 14 Hz, the longest 
time scale of variability - a few days at 10 12 Hz - is short enough 
to be consistent with the observations. Alternatively, variable op- 
tical depths within the emitting region could cause the density of 
SSC seed photons at frequencies ~ 10 12 Hz to fluctuate with- 
out changing the observed synchrotron radiation much, thereby 
causing the SSC X-rays to be more highly variable than the op- 
tical synchrotron radiation. 

The continuum spectrum from UV to near-IR wavelengths 
is well described by a power law with constant spectral index 
ttopt = 1.6, even as the flux varied by over an order of mag- 
nitude. As discussed above, this agrees with constant injection 
of relativistic electrons with a power-law energy distribution of 
slope -3.2 that is steepened to -4.2 at high energies by radiative 
energy losses. Although such steep slopes are not in conflict with 



particle acceleration models, neither are they predicted. Yet the 
constancy of this steep slope during the campaign implies that 
the physical parameters governing the particle acceleration pro- 
cess can be maintained over long periods of time and from one 
event in the jet to the next. This presents a challenge to theories 
for particle acceleration. 

The radio light curves from 230 to 1 GHz become pro- 
gressively smoother with decreasing frequency. As can be in- 
ferred from Fig. [10] there is also a lag in the timing of max- 
ima and minima in the light curves toward lower frequencies, 
as well as a delay of 100-150 days between optical variations 
and those at 5 GHz. This is readily explained as a consequence 
of opacity, with an outburst delayed at radio frequencies until 
the disturbance reaches the locati on in the jet where the optical 
depth is less than unity ( see, e.g., [Hughes, Aller. & AlleiLll985t 
Marscher & Gearl [T985). In addition, the radiative lifetime of 
electrons in the radio-emitting regions is longer than is the case 
farther upstream, owing to the lower magnetic field and photon 
energy density. This prolongs the duration of an outburst - and 
therefore the timing of the minima in the light curves - at lower 
frequencies. Because of the absence of optical and 230 GHz data 
near the global peak in 2006-07, we cannot determine whether 
there is a significant delay between the optical and 230 GHz light 
curves. However, the fact that at the end of the campaign period 
we see a brightening in the optical, but not in the mm light curve, 
suggests that the two emission regions are not co-spatial. 

8.2. Polarimetric behaviour 

The main feature of interest in the polarimetric data is the rota- 
tion of the fl-band and 43 GHz core EVPAs between MJD 54120 
and 54200 (see Fig. [9j. The optical polarization was several per- 
cent or less during the middle of this rotation, but it was as 
high as 23% near the beginning and end. Apparent rotations by 
~ 300°or more can occur simply from random walks caused by 
turbulent cells with random magnet i c field orientation passing 
through the emission region dJonesL [[988). 

In this case, how- 
ever, the apparent rotation has an ex tremely low probability t o 
be as smooth as we have observed dD'Arcangelo et afl, 12007). 
We therefore conclude that the rotation is intrinsic to the jet. The 
coincident rotation of the 43 GHz core and optical polarization 
vector suggests a common emission site at the two wavebands. 
[Note: the discrepancy between 43 GHz and optical EVPAs mea- 
sured at essentially the same time could be due to Faraday rota- 
tion, whic h is of the order of tens of degrees and migh t vary with 
location dZavala & Tavloill2004tlJorstad et al.Ll2007l) .l 

The observed behaviour is qualitat i vely similar to the mod - 
els suggested by IKikuchi et al.1 d!988l) : ISillanpaa etaTI d 19931) : 
IMarscher et alj d2008l) . where a shock wave or other compres- 
sive feature propagating down the jet traces a spiral path and cy- 
cles through the orientations of an underlying helical magnetic 
field. This manifests itself through rotation of the position angle 
of linear polarization as the feature moves outward. The polar- 
ization is low during the rotation because the symmetry of the 
toroidal component of the helical magnetic field produces a can- 
cellation of the net linear polarization when integrated over the 
entire structure of the feature. On the other hand, the cancella- 
tion is only partial, which implies that the feature does not extend 
over the entire cross-section of the jet. 

Several document ed events of opt i cal po sitio n angle rotation 
have b een reported by Sillanpaa et al. ( 1993) and Mar scher et al.l 
d2008h (both - BL Lac). IKikuchi et all dl988l) (OJ287), and 
lLarionov et all d2008l) (S5 0716+71). The time scale of rota- 
tion of the polarization vector in 3C 279 is much longer than in 
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any of these previous cases - almost two months as compared to 
~ 1 week. This difference can be explained by the substantially 
larger scale of the core of the 3C 279 jet as well as the greater 
distance from the central engine (as deter mined by the intrin - 
sic opening angle of the jet of 0.4°+0.2°; iJorstad et all [2005). 
The time required for the disturbance to reach and pass through 
the millimeter-wave core is therefore longer than in the other 
blazars. 

The first superluminal knot that emerged during our cam- 
paign (see Fig. [8]) was located at the centroid of the 43 GHz 
core on MJD 53888+55. (It is interesting that this coincided 
with the end of the interval of rapid fluctuations of the optical 
and, especially, X-ray flux.) The flux subsequently rose at all 
wavebands, reaching a peak ~ 130 days later. We note that the 
EVPA of the core at 43 GHz changed by ~ 130° as the flux 
rose, although the data are too sparse to determine whether this 
was a smooth rotation. The next knot coincided with the core on 
MJD 54063 ±40, ~ 50 days after the radio and X-ray flux peak 
but at a time of multiple, major optical flares. This was shortly 
before the start of the optical EVPA rotation featured in Fig. [9] 
The timing suggests that the optical flares and EVPA rotation in 
early 2007 took place in this knot as it moved downstream of the 
core. This contrasts with the case of BL Lac, in which the re- 
gion with helical ma gnetic field lies upstrea m of the millimeter- 
wave core. This led iMarscher et al] d2008l) to associate the he- 
lical field with the acceleration and collimation zone (ACZ) of 
the jet flow, downstream of which the flow is turbulent. For the 
polarization rotation to occur downstream of the 43 GHz core, 
either the core lies within the ACZ in 3C 279 or the helical field 
ca n persist beyond the A CZ. The latter possibility is favoured 
bv lGabuzda et alj ([2008) based on circular polarization observa- 
tions of blazars. The downward curvature of the EVPA vs. time 
curve of Fig.[9]implies that the pitch angle of the helix becomes 
smaller - i.e., the helix opens up - with distance from the core. 
The time scale for this to occur ~ 60 days, which corresponds to 
a distance of ~ 20 pc along the jet of 3C 279, given the kinemat- 
ics and angle to the line of sight of the jet derived in Sect. [4] 

9. Conclusions 

By following the evolution of the flux at radio, near-IR, opti- 
cal, UV, and X-ray frequencies and the linear polarization at ra- 
dio and optical bands intensely over a two-year period, we have 
uncovered patterns that reveal key aspects of the physics in the 
relativistic jet in 3C 279. The IR-optical-UV continuum spec- 
trum of the variable component follows a power law with a con- 
stant slope of -1.6, while that in the 2.4-10 keV X-ray band 
varies in slope from -1.1 to -1.6. This agrees with the expec- 
tations of an emission region into which electrons are steadily 
injected with a power-law energy distribution of slope -3.2 that 
is modified to a slope of -4.2 at high energies owing to radiative 
losses. The steepest X-ray spectrum occurs at a flux minimum. 
The least contrived explanation is that the X-ray emission at this 
time comes from a new component in an upstream section of the 
jet where the radiative losses from inverse Compton scattering of 
seed photons from the broad emission-line region are important. 
If this is the case, then a y-ray flare should precede the rising por- 
tion of a multi-waveband outburst, a prediction that can be tested 
with GLAST and AGILE along with intensive multi-waveband 
monitoring. 

During the decline of flux from the maximum in early 2007, 
we observe a rotation of the optical and 43 GHz core polarization 
vectors totaling ~ 300°. The smoot hness of the ro t ation leads 
us to conclude that, as in BL Lac (IMarscher et all 120081) and 



possibly other blazars, the jet contains a helical magnetic field. 
However, in contrast with BL Lac, the region of helical field in 
3C 279 extends ~ 20 pc past the 43 GHz core. Given the paucity 
of well-sampled optical polarization monitoring over periods of 
time longer than 1-2 weeks, such rotations may be the norm 
in blazars. If so, more extensive future polarization monitoring 
should uncover many more examples, allowing more general in- 
ferences to be drawn regarding the magnetic field structure in 
jets and its relationship to the dynamics of the flow. 
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